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We consider a nematic liquid crystal confined by two parallel planar interfaces, one being laterally 
homogeneous and the other provided by a substrate endowed with a periodic chemical stripe pattern. 
Based on continuum theory we analyze the influence of the lateral pattern on the liquid-crystalline 
Casimir force acting on the interfaces of the nematic cell at distance d due to the thermal fluctuations 
of the nematic director. For d much larger than the pattern periodicity, the influence of the patterned 
substrate can be described by a homogeneous, effective anchoring energy. By tuning this parameter 
we recover previous results for the liquid-crystalline Casimir force. For the general case, i.e., smaller 
separations, we present new numerical results. 
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I. INTRODUCTION 



Liquid crystals in general are sensitive over a wide 
spatial range to the anchoring conditions of confining 
interfaces. This holds also for lateral variations of an- 
choring conditions generated either by surface topo gra- 

$ffM MMMJIL U&M or Patterning H H [H 111 
[131 ll-l Ha Ha. Il7l Hq giving rise to numerous possible 
applications. In recent years, the influence of the struc- 
tured substrates on the properties of liquid crystals has 
been studied and it has been demonstrated that to a cer- 
tain extent such nontrivial geometries may optimize the 
performance of electro-optical liquid-crystalline devices. 
For instance, a four-domain twisted nematic liquid crys- 
tal display provides a wide viewing angle with no gray 
scale inversion )l9j and using multistable nematic liquid- 
crystal devices with micropatterned substrate alignments 
reduces the energy consumption 0, |(| 0, 0] . 

The influence of the anchoring on the liquid crystal 
order parameter translates into an effective interaction 
between the confining interfaces which may be provided 
either by true solid substrates or by an adjacent vapor 
phase where the former case can support permanent lat- 
eral structures. Here we consider liquid crystals deep in 
their nematic phase where the orientational order is de- 
scribed by a director field with long-ranged correlations. 
In case of any mismatch between the prescribed align- 
ment of the bulk and the substrates, the director struc- 
ture may not be uniform. In such a case, the free en- 
ergy of the system typically exhibits several, metastable, 
minima which upon a change in the parameters of the 
system - such as the film thickness, external, or internal 
forces - may turn into the glob al minimum resulting in a 
structural phase transition |20j . This makes the stability 
of the equilibrium configuration geometry dependent. In 
this context one should keep in mind that perturbative 



approaches may miss the occurrence of first-order struc- 
tural phase transition. We shall consider the case of frus- 
trated systems in which, however, the director structure 
remains uniform up to a critical thickness d c . Within 
such a regime the liquid-crystalline mediated effective 
interaction between the confining interfaces due to the 
thermal fluctuations of the director adds to background 
contributions due to structural forces arising from pres- 
mectic layerin g [2ll l22l |23| and enhanced ordering near 
the substrates |21| and due to dispersion forces |24j which 
exhibit only a weak temperature dependence. 

Using the continuum Frank free energy, we study 
the fluctuation-induced interaction - the so-called liquid- 
crystalline Casimir effect - between two parallel inter- 
faces where one is periodically patterned and the other 
one is homogeneous. We consider a periodicity in the 
local anchoring energy and model the liquid crystal- 
substrate interaction by the Rapini-Papoular surface 
free energy. We investigate the modification of the 
fluctuation-induced force compared with its behavior for 
substrates with uniform anchoring conditions |25| as a 
function of the pattern periodicity £ and the character- 
istic length of the pattern ( a (see Fig.^). 

Two model systems are considered. One consists of a 
substrate with a pattern characterized by homeotropic 
anchoring of alternating strengths facing a second sub- 
strate at a distance d which is characterized by a uni- 
formly strong homeotropic anchoring. For this system 
the mean director is constant at any separation. By 
changing the boundary condition via the change of the 
patterning ratio Ca/Ci the character of the force changes. 
Depending on whether the boundaries are effectively 
similar-nonsimilar or similar-similar the force is repul- 
sive or attractive, respectively. For certain values of Ca/C 
and of the reduced distance d/(, the liquid-crystalline 
Casimir force vanishes. In the second system the pattern 
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consists of alternating stripes of homeotropic and degen- 
erate planar anchoring while the upper substrate still ex- 
hibits strong homeotropic anchoring. In this case there 
is the possibility of texture formation |l6j. However, for 
separations smaller than a critical one the director struc- 
ture is indeed uniform [2£j . For those ranges of the model 
parameters for which the director is constant (see, c.f., 
Subsec.l - 



. 1111 Bfl . we find a nonmonotonic behavior for the 
fluctuation-induced force. 

There are several techniques that can be used to 
create periodic anchoring conditions such as photo- 
alignment lla. the selectively thiol-functionalized photo- 
orientation , and atomic force lithography [jjj . In the 
latter case patterning at the nano-meter scale is reached 
which allows one to investigate more efficiently the in- 
fluence of the patterned substrate on the liquid crystal. 
The results of our study, for example, might be helpful in 
designing thin patterned liquid crystalline films for which 
the fluctuation-induced force plays a role for the stability 
of the film. 

In Sec.|n] we describe the system and the formalism 
that we apply for calculating the fluctuation induced ef- 
fective interaction. In Subsec. llII Al we investigate the 
force in the presence of a pattern of alternating anchoring 
strengths and in Subsec. llII Bl with a pattern of compet- 
ing anchoring conditions. We present analytical results 
for patterns of small and large scales and numerical re- 
sults for patterns at intermediate scales. In Sec. II VI we 
summarize our results. 



II. SYSTEM AND THEORETICAL MODEL 

We consider a liquid crystal in a nematic phase and 
confined by two parallel planar interfaces at a distance d 
(see Fig.^l. The liquid crystal is described by the Frank 
free energy [2tJ 



F = l I d 3 x 



K^V ■ nf + K 2 (n ■ V x n) 2 

+K 3 (n x V x n) 2 l , (1) 



where n denotes the director of the liquid crystal, V is the 
nematic volume, and K\, K 2l and K% are the splay, the 
twist, and the bend elastic constants, respectively. The 
interaction between the liquid crystal and the substrate is 
modeled by the Rapini-Papoular surface free energy |2£| 
given by 



F s = -~ / d 2 xTT(x)(n-e) 2 , 

1 J A 



(2) 



where x = [x, y) denotes the lateral coordinates of Carte- 
sian coordinates r = (x, z) , A is the surface area, W is 
the anchoring energy per area, and e is the easy direc- 
tion, i.e., the preferred direction of the director at the 
substrate if W > 0. For W < the director prefers 
the direction perpendicular to e. In the following we re- 
strict the discussion to the case that e is perpendicular 
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FIG. 1: The geometry of the nematic cell with cross section 
A and volume V = Ad. The upper boundary is characterized 
by strong homeotropic anchoring. The lower substrate is pat- 
terned. The pattern consists of periodic stripes of anchoring 
energies per area W a and Wb with widths £ a and £(,, respec- 
tively, so that £ = £a + C&- The easy directions (see Eq. at 
both boundaries are normal to the interfaces and we consider 
such values of W a and Wb for which the thermal average of 
the director field no is homogeneous. 



to both substrates, which leads to homeotropic anchor- 
ing for W > and to degenerate (i.e., with no preferred 
azimuthal angle) planar anchoring for W < 0. 



On the lower substrate located at z — 0, we assume 
that the anchoring energy varies periodically along the 
x-direction. The pattern consists of alternating stripes of 
anchoring energies per area W a and Wb- The substrate 
remains translationally invariant in the y-direction. On 
the upper substrate located at z = d we assume uniformly 
strong homeotropic anchoring (see after, c.f., Eq. ©). 
In general, the local director field is given by 



n(x, z) = n (x, z) + 5(x, z) 



(3) 



where no(x, z) is the thermal average of the director and 
<5(x, z) is the fluctuating part with vanishing thermal av- 
erage (5) = 0. 



A. Mean-field behavior 

First we discuss the mean-field solution no(x, z) of 
the director field. In the case of homeotropic anchor- 
ing everywhere on the lower substrate, the uniform solu- 
tion no = (0,0,1) is the equilibrium configuration. On 
the other hand in case of planar anchoring everywhere 
on the lower substrate, the liquid crystal is subject to 
competing surface interactions at the top and the bot- 
tom. For this so-called hybrid cell - a cell with uniform 
homeotropic and uniform planar anchoring on each sub- 
strate - it has been shown that the substrate whose an- 
choring is stronger can impose a uniform director config- 
uration up to a critical separation between the plates [29j . 



3 



In the case of periodic pattern of homeotropic and pla- 
nar anchoring on the substrate fSubsec. llll 13 [I . the full 
phase diagram of the system within mean field theory 
and the structural phase transition between a uniform 
director configuration and a distorted one can be studied 
by means of numerical minimization of the free energy 
functional 0, |2(J . However, one can naively expect that 
in this system the tendency to form a uniform director 
field is enhanced with respect to the hybrid cell due to 
the presence of the interlaced homeotropic stripes. In the 
following we restrict the discussion to separations smaller 
than the critical separation d c for which the director con- 
figuration is uniform and focus on the fluctuations. 



for the lower substrate. We assume homogeneous anchor- 
ing on the upper boundary so that 

K= d M = \w j A d 2 x („(x, z = d))\ (7) 

Minimization of the Hamiltonian TL[v] — HtmikM + 
^smf M + ^surf M leads to the following boundary con- 
ditions: 

— Kd z v(x, z) + W a v(x, z) a(x) 

+ W b v{x, z) (1 - a{x)) = 0, z = 0, (8 a) 



Fluctuations of the director 



Kd z v(x., z) + Wv{x, z) = 0, z = d, 



(8 b) 



Next we consider fluctuations around the uniform di- 
rector no = (0,0,1). Since the director n(x, z) is a 
unit vector the fluctuations can be described by d = 

(6 x ,S y ,-l + y/l-Sl- SI) ^ (S x ,S y ,-6 2 j2 - 81/2), 

where S x (x, z) and S y (x, z) are the two independent 
components. According to Eq. and within the one- 
constant approximation, the bulk contribution to the sta- 
tistical weight exp (— /37ibuik[n]) for a director configu- 
ration n is given by H bu ik[n] = Hbu\k[S x ]+Hbu\k[S y ] with 



WbulkM = -j J d 3 x[Vu(x,z)Y 



(4) 



where = fc^T is the thermal energy and K is the 
effective elastic constant [30J . Eq. amounts to con- 
sidering Gaussian fluctuations, i.e., the Hamiltonian is 
quadratic in z/(x, z). This is expected to give a qualita- 
tively correct description of the system except near an 
incipient structural phase transition. 

As a local contribution the surface interaction is eval- 
uated at the interfaces z = 0, z = d. The lower substrate 
is characterized by the patterning function 



a(x)= J2 © 



k( 



c, 



2 



(5) 



where 0(x) is the Heaviside step function, £ is the peri- 
odicity, £a is the width of the stripe characterized by W a , 
and Cfe — C ~ Ca i s f ne width of the stripe characterized 
by Wb (Fig.[TJ. The stripes form sharp chemical steps 
between them. The function a{x) is one at the regions 
characterized by W a and zero elsewhere. Accordingly, for 
this model the surface interaction [Eq. (J2J] disregarding 
constant terms is given by 7Y sur f [n] = 7i SU rl [<^x]+W sur f [Sy] 
with 



n z s =>] = \[Waj d 2 x (i/(x, z = 0)) 2 a(x) 



+W b / d 2 x(^(x,z = 0)) (1 -a(x)) 

J A 



(6) 



where v is either 5 X or S y . After integration by parts 
in Eq. Q and using the boundary conditions given by 
Eqs. (|8 a|l and 18b|l. the Hamiltonian Ti.[f] reduces to 
Tl[v] = — /yd 3 ii/(x,z)V 2 c(x, z). In terms of the so- 
called extrapolation lengths A a and A;, 



A a (fc) — 



K 



a(b) 



(9) 



the boundary condition [Eq. I|8 all ] on the patterned 
substrate z — reads -4i(x, z) = — \bd z v{x, z) + 

^1 + Xb ~ Xa a(x)j f(x, z) = 0. Assuming strong 

homeotropic anchoring W — oo (A = 0) at the upper 
boundary, Eq. I|8 b() leads to the Dirichlet boundary con- 
dition .42 (x, z = d) = z/(x, z = d) = 0. 

As an aside, we note the relation between the present 
model and those for surface critical phenomena. Rescal- 
ing the fluctuating field by K/ksT, the Hamiltonian in 
Eq. JU can be written as H b [(p] = k B T J v d 3 x \ (Vip) 2 
and the surface interaction in Eq. J7J) is represented by 
Tis[f] = k B T f A d 2 x^ip 2 , where c = W/K is the inverse 
extrapolation length of the critical order parameter pro- 
file at a surface 31]. The limiting cases c = oo and c = 
correspond to Dirichlet and v. Neumann boundary con- 
ditions, respectively. The bulk Hamiltonian for a system 
close to the critical point also includes the terms ^ip 2 
and -^V? 4 °f which the former one vanishes at the criti- 
cal point. In this sense the present study corresponds to 
discussing, within the Gaussian approximation, a slab of 
a system at bulk criticality confined by planar surfaces 
one of them endowed with a pattern of the extrapolation 
length. 



C. The fluctuation-induced force 

We employ the path integral method introduced by Li 
and Kardar for calculating the partition function of the 
system [32, which amounts to integrate over all con- 
figurations of the fluctuating field weighted by the Boltz- 
mann factor and subject to the boundary conditions. We 
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impose the boundary conditions by inserting delta func- 
tions into the path integral. Thus the partition function 
Z of the field v reads 

Z = /2Mr) e- H M/^T) 

x n x 5(^i(x,z = 0))n x 5(^ 2 (x,^ = d)) (10) 

with the functional integral defined via a discretization 
on a lattice {r„} in the limit of a vanishing lattice con- 
poo di^(r„) 
— oo */2ir 

representation of the delta function, 



stant: / Vv(v) = JT„ f™ x 34]. Using the integral 



H x 8(Aa(xf) = j T>$ a exp(i J d 2 xt> a A a ^j , a = 1,2, 

(11) 



and performing the Gaussian integral over the field v, we 
obtain 



Z = M J II V ^ a e ' H °" ' 



(12) 



where ^0=1,2 are auxiliary fields defined at z = and 
z = d, respectively, W is a factor independent of d, and 
the effective interaction reads 

2 

H eS = I ^ x J d 2 x'* Q (x)M ct , /3 (x,x')* /3 (x'), 



a, 0=1 



where M is regarded as a matrix both with respect to 
the indices a, (3 and the coordinates x, x': 



\ h -\ 



M 11 (x,x') = {[l + ^a( 2 ;)][l , — 
M 12 (x,x') = (l- 



*o(ffi / )] + ^f^ [o(g) - a(z')R - A 2 3 2 }g(x -x', 2 - z') 

At, — A„ 



Z = 2'=0 



a(x') — Ab9 z 'JG(x — x', z — z') 
M 2 i(x, x') = (l + ^r*o(«) - A h c) z )G(x - x', z' - 2) 
Af 22 (x,x') = G(x-x',z-z') 
I 



=d.z'=0 



z'= d,z— 



(14) 



and the translational invariance along the y-direction, 
the matrix M in the lateral Fourier space, M(p, q) = 
/ / d 2 x d 2 x'M(x, x')e lp ' x e lqx , has the following form: 



where G(r,r') = i7r K B \J~r'\ 1S ^ ne bulk two-point correla- 
tion function in three dimensions defined by TpyrV 2 G(r— 
v') = -8{v-v'). 

Due to the symmetries in the xy-plane it is useful 
to switch to the Fourier transformed quantities. We 
note that if the patterning function a(x) = 1, i.e., for 
a homogeneous substrate (A a = A fc ), the matrix M 
is diagonal in the lateral Fourier space (p,q). How- 
ever, here the patterning function is piecewise either one 
or zero [Eq. JSJ] . Due to the periodicity of the pat- 
terning along the x-direction, i.e., a{x) — a{x + £), with the (2 x 2) matrices N m given by 



M(p,q) = (2ir) 2 5{p y + q y ) ^ N m {p xi p y ) 



m— — 00 



x 5 ( p x + q x + 



27rm 



(15) 



[(1 + ^t) 2 ~ X l d ^ G & z ~ z %=*>=o + ^ (1 + - A ^')G(p, z - z')\ z=d s/=0 

N = I 'I (16) 

(1 + _ Xbdz)G{pi z> z )\ ^ G(p, z - z%^ d 



and 



^2^a m [{l + ^*^)G(p, z-z') + G(p m , z - z')}\ z=z , =(j + </> m >^a m G(p, z - z')| 
x ~ K -a m G(p m ,z'-z)\ zl=dz=0 



=d,z'=0 



A„ 



with p = ijpl+p 2 ,, p m = ^j{ Px + 2Tim/C) 2 +p 2 y , the and 
two-point correlation function in lateral Fourier space / — A ^ 2 



G(p,z-z') = ^e-^- z \ (18) ' ' ' ' - 



a k a m - k G(p k ,0) (19) 
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where the prime at the summation sign indicates that 
in the sum the term k = is excluded. The patterning 
function a(x) is represented as a Fourier series a(x) = 
Efc-oo ^ 2mkxK with 

1 f L ' 2 , 1 

a/,. = — / dxa(x) expf— 2nikx/Q = — - sm.(Trk£ a /Q, 

L J_ L/2 TTfc 

where L is the extension of the system along the 
x— direction (see Fig.^J. We mention that the patterning 
function a(x) is coordinate dependent and the phase of 
the coefficients au depends on the choice of the position 
of the coordinate origin used for a(x), but as expected 
the final result for the force is independent of this choice. 
We have checked this numerically [c.f., Eq. (|25fl]. 

Each pair of the momenta (p, q) indicates one ele- 
ment of the matrix M which is a (2 x 2) matrix itself 
[Ea. I|15(l ]. Although M has an infinite number of ele- 
ments and is not diagonal it can be brought into a block 



diagonal form by an even number of row and column 
permutations [35l |3(j . We take the system to be peri- 
odically extended with period L = NC, along the cc-axis 
with N being an integer. This leads to momenta p x that 
are integer multiples of 2ir/L. For p x fixed the factor 
o(p x + q x + 2irm/C,) in Eq. I|15|l leads to nonvanishing 
matrix elements for all q x — ~p x ~ 2irm/( > with m G Z. 
This allows one to identify the block structure of the ma- 
trix M. For fixed j the momenta p x = 2irj/L + 2wl/( 
and q x — —2irj/L — 2nk/( t form the block Mj where 
l,k G Z and j — 1, ■•• ,N = L/C, so that there is no 
multiple counting of the momenta. One can read off 
the elements of the infinite-dimensional block matrices 
Mj(p yi q y ) from Eq. fT5|) : 

Mj,ki{p y ,q y ) = 2n5(p y + q y ) B k i (~J^>Py) ( 21 ) 
with the matrix B given by 



2vr 

B\p x = —3 



N ( Px -f) N_ x (p x ) N_ 2 (p x + f) 
Nx( Px -^) N ( Px ) N^( Px + f) 
N 2 (p x -f) N x {j> x ) N { Px + f) 



(22) 



so that the matrix element Bki reads 

[~ ,Pv ) = Nm=k - 1 [~ + —>Py J ■ ( 23 ) 

Note that a reindexation of the (p yi (7 y )-subspace is not 
necessary as M is diagonal with respect to it [Eq. 1(15)1]. 
Now the value of the path integral [Eq. I|12f) ] given by 

Z = N(det AfT 1/2 (24) 

can be calculated in Fourier space. The free energy 
— ksTlnZ of the system usually lends itself to a de- 
composition into bulk, surface, and finite-size contribu- 
tions 37] . The term — ksT In Af leads to the bulk free 
energy. The factor M introduced in Eq. I|12l) is given by 
TV = [det G _1 (r, r')] -1 / 2 where the determinant of the 
inverse of the two-point correlation function G(r — r') 
in an unperturbed nematic is calculated in the space 
actually occupied by the nematic, i.e., the volume V. 
Thus the result for the bulk free energy is given by 

Fbuik = k B TVf* m " (^g 2 ln ( KQ k ^f 2 ) where Q max 
is an ultraviolet momentum cutoff of the order of the in- 
verse size of the nematic molecules. The remaining part 
of the free energy F = A(F*^+F*=f)+5F(d) is given by 
F = ^|2IlndetM. Here the (i-independent terms F*=f 



and F*^f are the surface tensions associated with the in- 
terfaces at z — and z — d, respectively, and the finite- 
size contribution 5F(d) is the fluctuation-induced inter- 
action. In the present model there are no other finite-size 
contributions. Using the block structure of M, we ob- 
tain In det M = Influx det Mj) = E^i ln ( detM i) • 
Thus the fluctuation- induced force T = —ddF equals 
T = - ! ^Y,f=i TT ( M i 1 ddM j ). The final result for 
the force reads 




x Tr (B- 1 {p x , Py )d d B(p x ,p y )) , (25) 

where we have carried out the thermodynamic limit 
L — > oo so that the summation over j is replaced by 
So'" ^P x ' trace over the continuous momenta 

p y is also replaced by — J °° dp y and A — L 2 . Equa- 
tion l|25(l takes also into account that there are two in- 
dependent fluctuating fields v (i.e., 5 X and S y ) which 
lead to a doubling of the force. Note that in Eq. Ij25(l 
the trace is taken with respect to the remaining discrete 
indices k and I [Eq. J53J)]. This can be calculated nu- 
merically by truncating the matrix B k i at order /, i.e., 
k, I = -(1-1) /2, . . . , 0, . . . , (I-l)/2. The force T follows 
from extrapolating / — > oo. 



6 



III. RESULTS 

In the limit d/C ^> 1 the contributions from the ma- 
trices N m to det M decrease rapidly with increasing m 
[Eq. i|15fl ]. In this limiting case it is sufficient to consider 
only the contribution from No(p x ,p y ), in the sense that 
truncating the matrix B [Eq. at J > 1 leaves the 
integrand in Eq. 1)25(1 practically unchanged. Therefore it 
is instructive first to focus on this limiting case, which 
may correspond to a nano-patterned substrate facing a 
homogeneous substrate at a micrometer separation, and 
to investigate analytically the behavior of the force. 



h 




A. Pattern of different anchoring strengths 

Here we consider a pattern characterized by stripes 
of homeotropic anchoring so that both W a and W& are 
positive. In this case the fluctuations are suppressed at 
the substrates [Eq. ©]. Using Eas. l|53|l. (HI , and 
we find for the force 



H d > c) 



lim 



-k B TA 
2vr 2 



d P y 



27T/C 



dp x 



x Tr (Nq \p x ,py)d d N ( Px ,p y )) 



k B TA 
nd 3 



dx 



x+d/K!! exp (2x) 



x — d/\ a {{ 



where 



"Vff 



CaAfc + CbA a 



(26) 



(27) 



is introduced as an effective extrapolation length and 
x = pd is the rescaled momentum. Thus in the limit 
d^> £ the patterned substrate can be described by an ef- 
fective anchoring energy per area with the force found be- 
tween two homogeneous substrates where one substrate 
is characterized by strong anchoring and the other sub- 
strate is characterized by a finite anchoring, i.e., a finite 
extrapolation length A e ff |38|. 



If X a or Ab is zero, A e ff vanishes and the force is long- 
ranged and attractive: 



T[d » C, A off - 0) 



k B TA C(3) 
iird 3 



(28) 



where £(s) = X)*=i ^ S ^ s the Niemann zeta func- 
tion. This expression equals the one obtained for 
substrates both characterized by homogeneous strong 
anchoring [25|. This implies that in this limit the 
fluctuation-induced force is not affected by the stripes 
with weaker anchoring. This behavior resembles the one 
of the electrodynamic Casimir force between a flat and 
a rectangularly corrugated substrate |35l ] in the limiting 



FIG. 2: The fluctuation-induced force T as a function of the 
reduced separation d/\ e s in the case in which the patterned 
substrate can be described by an effective homeotropic an- 
choring (d 2> Ci see Eq. l26H h A crossover from repulsion to 
attraction occurs for a separation d/X c g ~ 1.19 followed by a 



minimum at 



,il ' k B TA/\-t 



(1.65,-0.003). 



case that the periodicity is much smaller than the am- 
plitude 2h of the corrugation and the mean separation 
H between the substrates. In this case the force equals 
the one between two flat substrates at a reduced mean 
separation H — h so that the force is not affected by the 
valleys of the corrugated substrate. 

At intermediate values of A e ff there is a crossover from 
attraction to repulsion (Fig. 0). At separations smaller 
than A* ~ 1.19 A ff, the boundaries effectively act as 
being homogeneous but dissimilar - one boundary char- 
acterized by strong anchoring and the other boundary 
characterized by a finite weak anchoring so that the force 
is repulsive • The asymptotic behavior of the force for 
dj Aeff <C 1 is given by 



T{C, < d < Aeff) 



3fc s TAC(3) 
167rd 3 



1 - 



81n2 d 
3C(3) A^ 



Thus in this regime the leading long-ranged repulsion 
term <^d~ 3 , corresponding to two homogeneous sub- 
strates characterized by infinitely strong and zero anchor- 
ing (Dirichlct- Neumann boundary conditions), is weak- 
ened. At separations larger than A* , the boundaries effec- 
tively act as being similar - one boundary characterized 
by infinitely strong and the other boundary character- 
ized by finite yet strong anchoring; therefore the force is 
attractive [3^| . The asymptotic behavior of the force for 
d/Aeff ^ 1 is given by (compare Eq. (J2HJ) 



^■(d> A c ff,C) » - 



k B TA((3) 
ind 3 



_ 3A c ff \ 
d J 



(30) 



This means that the long-ranged attraction ~d~ 3 , cor- 
responding to two homogeneous substrates characterized 
by infinitely strong anchoring, is reduced. In Fig.[3| the 
amplitude of the fluctuation- induced force F(d ^> £), as 



7 




1 x 1CT 5 - 



0.0 



-1 x 10~ 5 - 



-2 x 10~ 5 t 



"1 1 1 1 

A a /C = 5,A 6 /C=12,d/C = 10 




0.2 0.4 0.6 0. 

Ca/C 



FIG. 3: Amplitude of the fluctuation-induced force T relative 
to asymptotic behavior d~ A as a function of the reduced sepa- 
ration d/Aefi in the case in which the patterned substrate can 
be described by an effective homeotropic anchoring (d 2> C 
see Eq. 12611 ). At d/A e g = this amplitude attains the value 
3C(3)/(16tt) ~ 0.072 linearly [Eq. jf^J] and approaches the 
value -C(3)/(4tt) ~ -0.096 for d/X cS -» oo [Eq. lEiul l. 

given by the full expression in Eq. 1|26[) . is shown as a 
function of the reduced separation d/X e g. 



FIG. 4: The fluctuation-induced force J as a function of the 
patterning ratio Ca/C for a fixed reduced separation d/C (ful- 
filling the requirement d ^> C, see Fig. |2J and fixed reduced 
extrapolation lengths \ a (b)/C fEq. 1261 1. The crossover from 
repulsion to attraction indicates that upon increasing Ca/C 
the system transforms from the effective strong-weak to the 
effective strong-strong anchoring regime. For the given set 
of the parameters (Aa/C = 5, \b/C = 12, and d/C = 10), 
the force becomes attractive when more than 31% of the sub- 
strate consists of strong anchoring parts with W a > Wb, i.e., 
A a = K/W a < Xb = K/Wb- Note that here T is measured in 
units of C 3 instead of A^g as in Fig.|2] 



In Fig.|3] the force [Eq. (|26| is shown as a function of the 
patterning ratio Ca/C- Also in this case the above consid- 
erations provide an understanding of the crossover from 
repulsion to attraction upon changing the patterned sub- 
strate from the effectively weak to the effectively strong 
anchoring regime. 

In the opposite limiting case d/C, <C 1, the force is given 
by 

71(1 \ c Ja i^ijt ex P ( 2x ) + 1 

+ ^ [°° dx +*n ^ V( 31 ) 

C-/o ^exp(2,) + iy 

Here the geometrically weighted average is analogous to 
the results obtained within the proximity force approx- 
imation 01 • This scheme amounts to using the force 
density obtained from the homogeneous case and inte- 
grating over the local contributions of the force densities 
corresponding to the regions with anchoring energies per 
area W a and W&, respectively. The particular contribu- 
tions from the regions close to the chemical steps cause 
deviations from this approximation. Therefore this result 
holds only for a low number density of chemical steps, i.e., 
for d < C 

For intermediate values of d/C, we have calculated the 
force T numerically based on the complete expression 
given by Eq. I|25|l . In order to highlight the deviation 



of the fluctuation-induced force in the considered pat- 
terned system from the long-ranged behavior in the ho- 
mogeneous cases, we show in Figs.|S] and the rescaled 
decrease of the force relative to the decay ~ d~ 3 in the 
repulsion and attraction regime, respectively. 



B. Pattern of competing anchoring energies 

In the following we generalize the system to a pat- 
terned bottom substrate characterized not only by dif- 
ferent anchoring strengths but also different preferred 
molecular axes. We consider a pattern consisting of alter- 
nating stripes of either homeotropic or degenerate planar 
anchoring, i.e., positive values for W a and negative values 
for Wb, so that the liquid crystal is subject to compet- 
ing preferred orientations at the bottom substrate while 
the upper substrate still exhibits strong homeotropic an- 
choring. In this case fluctuations are suppressed at the 
stripes of homeotropic anchoring and enhanced at the 
stripes of planar anchoring [Eq. @]. For those ranges of 
the model parameters for which the frustrating effect of 
the interlaced planar anchoring does not modify the mean 
orientation of the director, the uniform configuration of 
the director is thcrmodynamically stable and the system 
can be described as in the previous subsection. However, 
if the planar anchoring becomes dominant, the uniform 
director configuration is destabilized. For this case, we 
determine the behavior of the fluctuation-induced force 
before the onset of the ensuing structural phase transi- 
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FIG. 5: Rescaled relative decrease of the fluctuation-induced 
force in the case in which the pattern is characterized by al- 
ternating anchoring strengths as a function of the reduced 
separation d/Xb compared to the repulsive force J- Icp = 
3fcsT^4^(3)/(167T(i 3 ) between two homogeneous boundaries 
characterized by infinitely strong and zero anchoring energy, 
respectively. For d < ( Eq. 1311 holds. Both integrals in 
Eq. 1311 happen to have the same form as in Eq. JUJ (which 
holds, however, for d 2> C) with A e ff replaced by \ a and Xb, 
respectively. According to the discussion in the second para- 
graph following Eq. 1261 , for the parameter values Q a = C,b 
and A a /Ai, = 0.4 chosen here these integrals are positive 
for d/\b < 0.48 for the contribution associated with X a and 
d/Xb < 1.19 for the contribution associated with Ab. Thus 
one expects that on the basis of Eq. 1311 J- reduces to J- rep as 
given by Eq. 12911 for d/£ — > 0. The solid curve corresponds 
to T given by Eq. 13111 and the other curves correspond to the 
full numerical results obtained from Eq. 1251 for the indicated 
values of (,/Xb- All curves appear to vanish for d/C, — > and 
thus confirm the above expectation. Moreover, by increas- 
ing the periodicity the difference between the solid curve and 
the numerical results vanishes and thus confirms Eq. I3H as 
a reliable approximation in the limit d/C, — » 0. 



tion. 

We consider again first the case d/C, 3> 1 which allows 
us to describe the force in terms of an effective anchor- 
ing energy and an effective extrapolation length. Within 
the same approximation as in the preceding subsection, 
i.e., truncating the matrix B at order 1 = 1, the force 
[Eq. (123] is given by 



T(d > C) 



where 



k B TA 
nd 3 



d.r 



x-d/\ ot{ 
x+d/\ a a 



exp (2x) + 1 



A,, 



(32) 



(33) 



and x = pd is the rescaled momentum 1 1 it . The form of 
this force is that of the force for the hybrid cell with 



FIG. 6: Rescaled relative decrease of the fluctuation-induced 
force in the case in which the pattern is characterized by al- 
ternating anchoring strengths as a function of the reduced 
separation d/Xb compared to the attractive force J-^tti = 
— ksTA^Ci) / (Aitd' 3 ') between two homogeneous substrates 
characterized by infinitely strong anchoring energies. For 
d ^ C Eq. 1261 holds. According to the discussion in the sec- 
ond paragraph following Eq. 1261 this integral is negative for 
d/Xb > 0.68 for the parameter values £ a = (5 and X a /Xb = 0.4 
chosen here. Thus one expects that T reduces to T^ttr given 
by Eq. for (/d -> 0. The solid curve corresponds to T 
given by Eq. 1261 and the other curves correspond to the full 
numerical results obtained from Eq. 1251 for the indicated val- 
ues of C/Xb- All curves appear to vanish for C,/d — > and thus 
confirm the above expectation. Moreover, by decreasing the 
periodicity C, the difference between the solid curve and the 
full numerical results vanishes and thus confirms Eq. 1261 as 
a reliable approximation in the limit (,/d — » 0. 



the bottom substrate described by an effective planar 
anchoring W cS = ( a W a /( + ( b W b /( < 0. Thus in this 
limit the patterning enters only via this expression for 
W cS . 

For £ d <C A e ff the long-ranged repulsion between 
the boundaries, characterized by strong and effectively 
weak anchoring, is enhanced: 



^(C < d < A cff ) 



3k B TA((3) 
16rf 



1 



81n2 d 
3C(3) A^ 



(34) 

Upon approaching the critical separation d c = A e s the 
force increases and, within the Gaussian approximation, 
diverges logarithmically at d c : 



Hd/d c 



1) 



ik B TA 
^7 



In 1- 



■ In 2 



(35) 

This behavior follows from the fact that for d = A e fj the 
denominator of the integrand in Eq. I|32|) vanishes CIS 7j 
for x —> 0. Thus the pretransitional behavior (d — > d c ) 
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FIG. 7: The fluctuation-induced force as a function of the re- 
duced separation d/Xb in the case in which the pattern is char- 
acterized by competing anchoring energies. The solid curve 
shows the result for small periodicity (£ <C d) where the pat- 
tern can be described by an effective anchoring energy per 
area W e g = ( a Wa/( + ( b W b /( < [Eq. fHJ]. The dashed 
curve shows the result for large periodicity 3> d) for which 
the proximity force approximation is valid [Eq. 136H 1. The 
data points (+) show our full numerical results [Eq. for 
£/A[, = 1.5. This illustrates that as expected T approaches 
the proximity force solution (dashed curve) for d/C, —> and 
the effective anchoring solution (solid curve) for d/C, — > oo. 



of the force is related to the sigularity of the soft mode 
(x -> 0). 

The logarithmic divergence is a characteristic feature 
of the perturbative method [Eq. if applied to a system 
subject to structural phase transitions |42l l43l 144. l45|. 
Here the planar anchoring destabilizes the uniform struc- 
ture governed by the strong homeotropic anchoring at the 
upper boundary. Upon increasing the separation d the 
influence of the upper boundary decreases and the desta- 
bilizing effect of the substrate characterized by the pla- 
nar anchoring increases so that upon approaching d = d c 
the uniform structure becomes unstable |2^]. Beyond 
the critical separation, the director structure is no longer 
uniform and the Gaussian Hamiltonian given by Eq. Q 
does no longer describe the system. For d> d c one must 
consider the fluctuations around a nonuniform configura- 
tion no(r). It is possible and even likely that the actual 
critical film thickness for the structural phase transition 
is not given by d c = A c ff as obtained in the present per- 
turbative approach. 

In the regime of the uniform director configuration, 
for separations comparable to the periodicity we have 
calculated the force numerically. As for Eq. (|31|l , in the 
limit d/( <C 1 the force is given by the geometrically 



weighted average of the local force densities: 



* d3 \<J° H£«TWti 



d/\ a 

Cj o UA ^exp(2 2 ;) + l / 



da;- 



(36) 



In Fig.[7|the behavior of the fluctuation-induced force is 
shown for the three different regimes. 



IV. SUMMARY AND CONCLUSION 

Based on the Frank [Eq. 0J] and Rapini-Papoular 
[Eq. (J2J] expression for the bulk and surface free energy, 
respectively, we have studied by field theoretical tech- 
niques the contribution of Gaussian fluctuations of the 
orientational order to the effective force acting on the 
planar and parallel boundaries of a nematic film of thick- 
ness d. The upper boundary exhibits strong and ho- 
mogeneous homeotropic anchoring whereas the bottom 
boundary is characterized by a one-dimensional, steplike 
periodic modulation with period £ of either the strength 
or the orientation of the anchoring which is homeotropic 
or planar (see Fig. . The system parameters are chosen 
such that the thermal average of the nematic order con- 
figuration is spatially homogeneous throughout the film. 
Our main results are the following: 

(1) For d 3> C) the fluctuation-induced force is propor- 
tional to d~ 3 times a scaling function which attains the 
same form as in the case that the patterned bottom sub- 
strate exhibits a homogeneous anchoring energy per area 
W e f[ = CaW a /C + (bWb/( with the different anchoring 
energies per area W a (h) — K/X a ^) weighted according to 
their lateral geometric contribution (Eqs. (|26|l and i|27|) 
as well as ll-S'-'f and H33|) ). If a change of the pattern 
of the anchoring strength causes the corresponding ef- 
fective boundary condition on the patterned substrate to 
change from homeotropic weak to homeotropic strong an- 
choring, the force exhibits a crossover from repulsion to 
attraction (see Figs. I2HH1 and^J. For a pattern of com- 
peting anchoring energies in the case that the stripes with 
planar anchoring are dominant relative to the stripes of 
homeotropic anchoring, the fluctuation-induced force is 
repulsive and nonmonotonic. 

(2) For d <^ C, the proximity approximation is valid, 
i.e., the force is given by the geometrically weighted av- 
erage of the local force densities obtained from the ho- 
mogeneous substrates (Eqs. (|31fl and l|36|l). 

(3) For intermediate separations, we have examined 
the behavior of the force numerically. For a pattern of 
homeotropic anchoring with alternating strengths, the 
rescaled relative decrease of the fluctuation-induced force 
as a function of the reduced separation compared to the 
long-ranged force is summarized in Figs .[S] and for sev- 
eral values of the reduced periodicity. We note that by 
increasing the reduced periodicity the system approaches 
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the limit of the proximity force approximation given by 
Eq. IpH)) (solid line in Fig.0, while by decreasing the 
reduced periodicity the system approaches the limit of 
the effective anchoring approximation given by Eq. l|2rj|) 
(solid line in Fig.|HJ). 

(4) For a pattern of competing anchoring energies, the 
behavior of the force is summarized in Fig.[7| Also in 
this case the full numerical solution interpolates between 
the asymptotic formulae given for small d [Eq. I|36|l ] and 
large d [Eq. Q ], 

Besides the force due to fluctuations of the director 
in a nematic film, there are additional forces acting be- 
tween interfaces of the film such as the well-known van 
der Waals dispersion force [24[ so that the total force is 
the sum of the fluctuation-induced and the dispersion 
forces. In the case that the mean director is inhomoge- 
neous so-called structural forces appear in addition which 
are not present in the homogeneous case considered here. 
For a film geometry of thickness d and area A, the disper- 
sion force decays as —AH/ (Gird 3 ) where H is the so-called 
Hamaker constant. As implied by Eas. (|26|l and Q32JI. 



the magnitude of the fluctuation-induced force scales as 
ksTA/d 3 . Thus, apart from a numerical prefactor of or- 
der one, the overall ratio between the dispersion force and 
the fluctation-induced force is given by H/ksT. Since liq- 
uid crystals and glass substrates have typically compara- 
ble indices of refraction, the Hamaker constant is of the 
order 10 -21 J ^(| which is the same order of magnitude as 
ksT at room temperature. Therefore typically the dis- 
persion and fluctuation induced force in liquid crystals 
are of the same order of magnitude and thus compara- 
ble. We mention that for a more adequate estimate of 
the Hamaker constant the anisotropy of the liquid crys- 
tal has to be taken into account [43 • For larger thick- 
nesses the dis pers ion force decays faster (~ 1/e? 4 ) due 
to retardation [43 , so that for sufficiently thick films the 
fluctuation-induced force dominates. 

Finally we note that while the direct measurement of 
the fluctuation-induced force in liquid crystals has not 
yet been accomplished, these forces affect in a character- 
istic way the pattern formation of thin liquid-crystalline 
dewetting films [H WL . 
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